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ABSTRACT 


The  pre-shot  predictions  and  evaluations  of  ground 
motion,  stemming  and  containment  are  presented  and  compared 
with  the  available  observed  data. 

The  analysis  of  seismic  data  included  corrections  for 
frequency  response  of  the  instruments,  derivation  of  non- 
recorded  motions  by  differentiation  and  integration,  and 
filtering  of  the  seismic  signals  to  derive  amplitude-frequency 
relationships.  Seismic  data  were  available  from  stations 
located  between  2  to  32  kilometers  from  surface  zero.  Ob¬ 
served  ground  motions  were  lower  than  expected  from  a 
coupled  (fully  tamped)  explosion  of  equal  size.  The  peak 
accelerations  and  velocities  displayed  high  rates  of  at¬ 
tenuation  compared  to  Salmon.  Relatively  small  azimuthal 
variations  were  observed.  Peak  motions  were  associated  with 
frequencies  between  3  and  80  hz .  A  noticeable  decrease  in 
frequency  with  increasing  distance  was  observed.  The  pur¬ 
pose  of  Sterling  was  to  test  decoupling  theory  which  pre¬ 
dicts  that  seismic  signals  can  be  reduced  by  detonating  a 
nuclear  device  in  a  large  cavity.  The  observed  decoupling 
effects,  were  strongly  frequency-dependent.  At  frequencies 
between  2  and  5  hz ,  particle  velocities  were  reduced  by  a 
factor  of  about  200.  Outside  this  range,  observed  decoupling 
factors  decreased  exponentially  with  frequency.  The  observed 
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frequency  dependence  of  the  decoupling  factor  was  found  to 
agree  with  decoupling  theory. 

The  prediction  of  containment  and  evaluation  of  the 
stemming  plan  were  substantiated  by  the  test  results.  No 
other  post -shot  data  were  available  for 


comparison . 


CHAPTER  1 


INTRODUCTION 


1 . 1  BACKGROUND 

Environmental  Research  Corporation,  under  contract  to 
the  Nevada  Operations  Office  of  the  U.  S.  Atomic  Energy 
Commission,  is  responsible  to  the  Safety  Evaluation  Division 
of  that  office  for  providing  safety  evaluations  of  under¬ 
ground  nuclear  detonations.  Specific  tasks  for  which  this 
organization  is  responsible  include  predictions  of  close- 
in  phenomenology  as  it  relates  to  containment,  overall  con¬ 
tainment  evaluations,  predictions  of  ground  motion,  and  re¬ 
duction  and  analysis  of  seismic  data. 

The  Sterling  event  took  place  on  December  3,  1966  at 
06:15  hours  CST  in  the  Tatum  Salt  Dome  in  Lamar  County, 
Mississippi.  The  nuclear  device  was  designed  to  have  a 
yield  of  0.36  kilotons  and  a  maximum  probable  yield  of  0.44 
kilotons .  The  actual  yield  was  estimated  to  be  0.35  kt 
(Reference  1.1).  It  was  located  at  a  depth  of  828  meters  in 
the  center  of  the  cavity  created  by  the  Salmon  event.  De¬ 
tailed  descriptions  of  the  Salmon  cavity  can  be  found  in 
References  1.2  and  2.10. 

The  purpose  of  the  Sterling  experiment  was  to  test 
decoupling  theory  which  states  that  seismic  signals  can  be 
reduced  by  detonation  of  the  device  in  a  cavity.  The 
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decoupling  factor  is  defined  as  the  ratio  of  the  ground 
motion  amplitudes  (particle  acceleration,  displacement,  and 
velocity)  from  a  tamped  (fully  coupled)  event  to  the  ampli¬ 
tudes  from  a  decoupled  event  of  equivalent  yield.  The 
amount  of  decoupling  is  dependent  on  the  ratio  of  the  yield 
to  the  cavity  volume,  and  to  some  extent  on  the  frequency 
of  the  generated  seismic  signals.  The  degree  of  decoupling 
expected  from  the  Sterling  event  was  calculated  by  the 
Lawrence  Radiation  Laboratory  to  be  between  factors  of  20 
and  160,  depending  on  the  amount  and  the  effects  of  the 
Salmon-caused  alteration  of  the  material  forming  the  cavity 
wall  (Reference  1.2). 

1.2  OBJECTIVES 

The  objectives  of  ERC i s  participation  in  Project  Sterling 
were  to  provide  the  Safety  Evaluation  Division  of  the  USAEC/NVOO 
with : 

1.  Pre- shot  predictions  of 

a.  Close-in  phenomenology  affecting  containment. 

b.  Peak  ground  motions. 

2.  Post-shot 

a.  Reduction  and  analysis  of  seismic  data. 

b.  Comparisons  of  predictions  and  observed 
results . 


-2- 


1.3  THEORY 


The  evaluation  of  containment  was  based  on  procedures 
developed  as  part  of  a  Long-Range  Program  directed  to  the 
continuing  improvement  of  confidence  in  predictions.  Cur¬ 
rently,  greater  dependence  is  placed  on  empirical  relation¬ 
ships,  and  it  is  probable  that  experience  will  always  be  a 
major  factor  in  evaluation.  However,  a  large  part  of  the 
improvement  will  come  from  greater  knowledge  of  the  pressure¬ 
time  history  and  temper ature-time  history  of  underground 
explosions.  This  is  being  pursued,  and  some  of  the  results 
were  employed  in  the  development  of  predictions  and  the 
evaluation  for  Sterling. 

The  predictions  of  ground  motion  were  based,  to  a 
large  degree,  on  experience  gained  from  the  Salmon  event 
which  was  detonated  at  the  same  location.  Except  for  the 
cavity  and  the  alteration  of  the  rock  adjacent  to  the 
cavity,  the  geologic  environment  affecting  transmission  of 
seismic  energy  from  the  Sterling  event  was  identical  with 
that  of  the  Salmon  event.  Therefore,  the  predicted  attenua¬ 
tion  of  amplitude  with  distance  was  taken  from  the  regression 
lines  through  the  data  from  Salmon. 

The  Sterling  yield  was  accounted  for  by  use  of  cube- 
root  scaling  relationships  to  scale  from  Salmon  (5  kt)  to 
Sterling,  .36  kt .  This  was  considered  appropriate  because 
the  predictions  for  Salmon  (which  agreed  with  the  data) 
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were  made  using  the  same  procedure  to  scale  from  Gnome, 

3  kt,  to  Salmon.  The  remaining  consideration  was  the 
amount  of  decoupling  which  might  occur  as  a  result  of  the 
Sterling  event  being  detonated  in  the  Salmon-created  cavity. 
Scaling  from  Salmon  provides  predictions  which  would  be 
applicable  only  if  Sterling  were  fully  coupled.  Because 
the  amount  of  decoupling  was  the  purpose  of  experiment  it 
was  not  possible  to  know  in  advance  how  much  decoupling 
might  occur .  In  order  to  make  predictions  which  would  be 
more  realistic,  the  fully  coupled  predictions  were  divided 
by  a  decoupling  factor.  The  Sterling  Technical  Concept 
(Reference  1.2)  stated  that  the  decoupling  factor  was  esti¬ 
mated  to  be  between  20  and  160,  depending  on  the  alteration 
of  the  cavity  wall  caused  by  Salmon.  To  provide  Safety 
Evaluation  Division,  NVOO,  with  ground  motion  estimates 
which  would  account  for  decoupling  effects  and  at  the  same 
time  would  be  conservative  for  safety  purposes,  the  small¬ 
est  decoupling  factor  (20)  given  in  the  Technical  Concept 
was  used.  Peak  accelerations,  displacements,  and  velocities 
predicted  for  the  tamped  (fully  coupled)  situation  were 
divided  by  this  decoupling  factor  to  provide  the  predicted 
peak  motions  for  the  decoupled  Sterling  event.  No  attempt 
was  made  to  consider  the  frequency  dependence  of  the  de¬ 
coupling  factor  in  the  pre-shot  predictions. 
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CHAPTER  2 


PROCEDURE 

2.1  GEOLOGIC  ENVIRONMENT 

2.1.1  Location  and  Topography 

The  site  of  the  Sterling  event  was  within  the  Tatum 
Salt  Dome,  Lamar  County,  Mississippi.  The  Tatum  surface 
area  is  about  75  to  100  meters  above  sea  level.  The  sur¬ 
face  is  moderately  dissected  with  narrow,  flat-topped  ridges 
rising  about  30  meters  above  intervening  valleys.  The  hills 
are  well  drained,  but  usually  the  bottomlands  are  wet  during 
most  of  the  year.  The  principal  streams  in  the  area  are 
Half  Moon  and  Grantham  Creeks  which  flow  into  Lower  Little 
Creek  about  one  or  two  kilometers  north  of  the  dome  (Refer¬ 
ence  2.1).  The  area  immediately  over  the  dome  is  a  topo¬ 
graphic  low. 

2.1.2  Regional  Geology 

Formations  outcropping  in  the  state  of  Mississippi  are 
of  upper  Cretaceous  and  Tertiary  age.  The  older  formations, 
exposed  along  the  northern  edge  of  the  state,  are  overlain 
to  the  south-southwest  by  progressively  younger  beds  which 
outcrop  in  roughly  parallel  bands.  In  general,  the  forma¬ 
tions  thicken  to  the  south-southwest.  Younger  beds  dip 
about  1.9  meters  per  kilometer  to  the  south,  whereas  dips 
on  the  base  of  the  older,  deeper  Eocene  beds  are  about  6.6 
meters  per  kilometer  southward  (Reference  2.1). 
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Tatum  dome  is  a  large  buried  salt  stock  that  has  pierced 
sedimentary  deposits  of  early  Cretaceous  to  Oligocene  (middle 
Tertiary)  age.  Sediments  of  Miocene  age  surround  and  overlie 
the  dome,  but  are  thinner  over  the  dome  than  in  the  surround¬ 
ing  area.  Rocks  of  Eocene  age  are  absent  on  the  crest  of  the 
dome  (Reference  2.2).  Regional  strikes  and  dips  are  locally 
disturbed  around  the  dome. 

2.1.3  Local  Stratigraphy  and  Lithology 

Stratigraphic  relations  over  and  around  the  dome  are 
shown  in  the  attached  cross  section  (Figure  2.1),  which  ex¬ 
tends  southwest  through  the  dome.  Figure  2.2  shows  the  loca¬ 
tion  of  this  section. 

Figure  2.3  is  a  generalized  columnar  section  at  Station 
1A,  the  location  of  the  Sterling  event.  Abbreviated  litho¬ 
logic  descriptions  of  the  various  rock  types  in  the  dome 
area  are  given  on  this  figure.  The  sediments  over  the  dome 
are  mainly  sands,  clays  and  silts  in  varying  mixtures.  The 
true  caprock  which  underlies  the  Catahoula  Sandstone  is  com¬ 
posed  of  two  dominant  lithologic  units.  The  upper  unit  is 
fine-  to  medium-grained,  crystalline  limestone.  The  lime¬ 
stone  contains  zones  of  high  permeability  in  which  varying 
quantities  of  drilling  fluid  were  reported  lost  while  drill¬ 
ing.  It  also  contains  some  lenses  of  medium-grained,  cal¬ 
careous,  loose  sand.  This  sand  may  be  either  contained  as 
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U.S.  DEPARTMENT  OF  INTERIOR  TECHNICAL  LETTER 

GEOLOGICAL  SURVEY  DRIBBLE- 29 


Adopted  by  U.S.  Geolooicol  Survey  from 
“Composite  Site  Mop',’  Project  Oribble, 
Oecember  4,1362,  by  Holmes  6  Norver,  Inc. 


Figure  2.2  Map  of  Tatum  Dome  Area  Showing  All  Holes  Drilled  on 
Dome  (Note:  A~Af  is  location  of  cross  section,  Figure  2.1) 


-3- 


-9- 


Figure  2.3 


lenses  within  the  limestone  or  as  inter finger ing  lenses  of 
the  Catahoula  Sandstone  with  the  caprock  limestone  (Refer¬ 
ence  2.3).  In  some  places,  the  limestone  has  been  leached, 
forming  channels  and  vugs.  The  lower  6  meters  of  the  lime¬ 
stone  (in  WP-4)  consist  of  vuggy  brecciated  black  calcitic 
limestone . 

The  limestone  at  its  base  grades  into  several  feet  of 
gypsum  which  in  turn  grades  down  into  the  anyhydrite  which 
composes  the  lower  dominant  lithologic  unit  of  the  caprock. 

The  anhydrite  is  hard  and  finely  crystalline.  In  WP-4,  the 
upper  part  of  the  anhydrite  has  been  highly  fractured,  and 
gypsum  veins  line  the  slickensi ded  surfaces.  Its  basal  con¬ 
tact  with  the  salt  is  sharp.  The  salt  varies  from  anhydrite¬ 
bearing  halite  to  transparent,  very  coarse  crystalline,  al¬ 
most  pure  halite.  The  salt  is  banded  with  nearly  vertical 
alternating  thin  bands  of  anhydrite-bearing  halite  and  al¬ 
most  pure  halite  (Reference  2.4). 

Tabulations  of  Pre-Salmon  drill  hole  information  and 
lithologic  logs  are  shown  in  Tables  3.3.1  and  3.3.2  of  Refer¬ 
ence  2.5.  Pre-Salmon  holes  are  shown  on  Figure  2.2.  Post- 
Salmon  holes  are  shown  on  Figure  2.4  (taken  directly  from 
Reference  2.10)  and  are  discussed  in  Section  2.2. 

2.1.4  Physical  Properties 

Physical  properties  for  the  salt  and  surrounding  sedi- 
nents  are  compiled  in  Tables  2.1  and  2.2.  Changes  in  physical 
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Unmarked  means  are 


properties  in  the  immediate  area  of  the  existing  cavity  are 
discussed  at  length  in  Reference  2.10. 

2.1.5  Structure 

The  enclosed  cross  section  (Figure  2.1),  extending 
southwest  through  HT-1  and  HT-2,  shows  the  structural  rela¬ 
tions  across  the  dome. 

Faulting  of  the  caprock  is  indicated  in  the  vicinity  of 
E-7  (See  Figure  2.1).  There,  a  41  meter  section  of  brecciated 
sandstone,  calcite  conglomerate  and  black  chert  gravel  was 
found  where  the  anhydrite  would  normally  be.  As  previously 
mentioned,  numerous  fractures  are  present  in  the  anhydrite. 

For  a  discussion  of  fractures  due  to  the  Salmon  event,  see 
Reference  2.10  and  Section  2.2  of  this  report. 

Although  little  mention  of  jointing  is  made  specifi¬ 
cally  in  the  literature,  it  seems  reasonable  to  expect  its 
moderate  development  in  the  more  brittle  lithologic  types 
including  sandstones,  limestones,  anhydrites  and  gypsum. 

The  salt  within  the  dome  is  essentially  free  from  jointing 
and  fracturing,  with  the  exception  of  Salmon-induced  fractur¬ 
ing  as  mentioned  above. 

2.1.6  Hydrology 

A  complete  review  of  the  hydrology  of  the  dome  area  is 
contained  in  Reference  2.2  (Pre-Salmon)  and  Reference  2.7 
(Post -Salmon ) . 
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Post-shot  data  indicate  that  the  Salmon  event  was  es¬ 


sentially  contained  within  the  salt  matrix  of  the  (Tatum) 
dome  and  no  explosion-produced  radioactivity  had  been  de¬ 
tected  in  the  aquifers  (Reference  2.7). 

2.2  PHENOMENOLOGY  AND  CONTAINMENT  PREDICTIONS 

2.2.1  Cavity  Radius 

The  device  for  the  Sterling  event  was  placed  in  the 
center  of  the  approximately  spherical  Salmon  cavity.  The 
horizontal  radius  of  the  cavity  was  about  16.5  meters 
(Reference  2.8). 

It  was  assumed  that  the  total  energy  for  this  event 
would  be  suddenly  distributed  uniformly  throughout  the  volume 
of  the  chamber.  Therefore,  the  resulting  initial  equili¬ 
brium  pressure  on  the  chamber  walls  could  be  estimated  on  the 
basis  of  the  work  done  in  an  adiabatic  process.  Thus 

p  =  .(Q/.~1.)w.  =  189  bars 
V 

where  a  =  1.2  for  air,  W  =  .44  kt  and  V  =  1.96x10^  cm3. 

The  total  heat  in  the  explosion,  based  on  a  yield  of 
.5  kt ,  is  4.4x10^  calories.  Neglecting  the  very  small  con¬ 
tribution  from  the  metal  of  the  device  and  any  steam  from 
contained  water,  the  energy  density  was  computed  as  follows: 

Yield  =  _ 4.-4.x\°.11 _  =  1.74xl04  cal/gm 

VxcYair  1 .96xl010xl  .293xl0~3 
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This  corresponded  to  a  temperature  of  22000°K  (see  Appendix 
A,  Table  A1 ,  Internal  Energy  Table).  An  indication  of  the 
rate  of  temperature  changes  for  a  .5  kt  explosion  in  a  14.2 
meter  radius  cavity  was  obtained  from  Appendix  D  of  Refer¬ 
ence  2.9.  Conditions  for  this  calculation  were  considered 
to  be  close  enough  to  those  for  Sterling  to  be  indicative. 

The  vaporization  of  the  salt  would  cause  a  rapid  drop  in 
temperature.  Cooling  would  occur  at  a  rate,  such  that  the 
vaporized  salt  would  condense  in  about  60  seconds,  based 
on  the  vaporization  temperature  of  2500°K  stated.  From 
that  time,  continued  melting  would  reduce  the  temperature 
to  the  melting  point  of  1059°K.  Solidification  would  begin 
at  about  2  hours  (Figure  A1 ,  Appendix  A)  and  from  that  point 
cooling  would  proceed  slowly. 

Computation  shows  that  281  cal/gm  is  required  to  raise 
salt  to  the  melting  point  and  to  melt  it.  Deducting  the  in¬ 
ternal  energy  in  the  cavity  at  the  melting  point  and  assuming 
that  the  remaining  heat  would  be  used  in  melting,  it  was 
estimated  that  the  thickness  of  salt  removed  would  be  20  cm. 

At  initial  high  pressures,  some  crushing  and  spalling 
of  the  salt  in  the  cavity  walls  was  anticipated.  This,  added 
to  melting,  was  predicted  to  result  in  a  total  increase  in 
cavity  radius  of  about  30  cm. 


2.2.2  Relationship  of  Internal  and  Over¬ 
burden  Pressure 

The  maximum  theoretical  pressure  of  189  bars  slightly 
exceeded  the  overburden  pressure  of  182  bars.  However,  it 
was  felt  that  if  this  theoretical  maximum  were  achieved,  it 
would  exceed  overburden  pressure  only  for  a  very  short  time. 
The  temperature  drop  to  reduce  the  pressure  from  189  bars 
to  182  would  be  about  1600°K.  Figure  A2,  Appendix  A,  indi¬ 
cated  that  such  a  reduction  in  temperature  would  occur  within 
milliseconds.  As  further  check,  a  calculation  was  made  of 
the  pressure  following  condensation  of  the  salt.  At  this 
time  all  the  metal  from  the  canister  would  have  condensed 
as  well,  since  the  vaporization  temperature  is  higher  for 
steel  than  salt.  Using  the  energy  density  for  air  at  2500°K, 
namely,  554  cal/gm,  the  pressure  was  computed  to  be  6  bars. 

As  indicated  previously,  this  condition  would  be  reached  in 
about  60  seconds. 

From  the  above,  it  was  apparent  that  the  pressure  in 
the  cavity  would  be  reduced  below  overburden  almost  instan¬ 
taneously  and  to  a  very  low  pressure  in  about  one  minute. 
Furthermore,  the  largest  theoretical  pressure  v;as  pre¬ 
dicted  without  consideration  for  any  energy  loss  by  vapori¬ 
zation  and  melting  of  halite.  With  conservative  considera¬ 
tion  of  energy  loss,  a  realistic  estimate  of  the  peak  pressure 
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indicated  that  it  would  probably  not  exceed  the  overburden 
pressure . 

2.2.3  Cracking  Radius 

Observed  fractures  from  the  Salmon  event  were  found  to 
extend  120  meters  from  the  weapon  point  (Reference  2.10). 

The  fracture  radius  of  Sterling  was  expected  to  be  much  less, 
and  extension  of  Salmon  fractures  by  the  Sterling  event  was 
not  anticipated. 

2.2.4  Radius  of  Radiation  Injection  into  Cracks 

Radioactivity  was  found  at  37  meters  from  the  Salmon 

W.P.  (Reference  2.10).  Extent  of  Sterling  radioactivity  was 
expected  to  be  much  less. 

2.2.5  Surface  Spalling 

The  depth  of  spall  for  a  tamped  event  of  the  yield  of 
Sterling  would  be  considerably  less  than  the  12  meters 
measured  at  Salmon  (Reference  2.11)  because  of  the  reduced 
motion  ax  ground  zero.  With  a  decoupled  event  the  motions, 
and  consequently  the  depth  of  spall,  would  be  reduced  further. 
Therefore,  spall  depth  was  considered  insignificant  to  con¬ 
tainment  for  Sterling. 

2.2.6  Chimney  Height 

Neither  Gnome  nor  Salmon  produced  a  chimney.  Therefore, 
no  chimney  was  expected  from  Sterling. 
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2.2.7  Analysis  of  Stemming 

Stemming  for  the  Sterling  emplacement  hole  is  shown  in 
Figure  2.5.  The  hole  was  filled  with  normal  sand  stemming, 
reinforced  at  the  bottom  by  a  157  meter  grout  plug. 

Because  of  the  length  of  this  plug  in  relation  to  the 
small  magnitude  of  close-in  deformation  effects,  it  seemed 
certain  that  the  plug’s  integrity  would  be  essentially  un¬ 
affected  by  the  blast.  As  a  safeguard  against  the  unlikely 
possibility  of  minor  seepage,  a  1.5  meter  thick  charcoal 
layer  was  placed  immediately  above  the  plug  to  diminish  any 
escaping  radioactivity.  A  similar  charcoal  layer  was  placed 
near  the  top  of  the  stemming  with  1.5  meters  of  impermeable 
polymer  above  it. 

In  light  of  the  conservative  features  of  its  design, 
the  stemming  was  considered  fully  adequate. 

The  two  post-shot  Salmon  drill  holes  were  completely 
grouted,  eliminating  the  hazard  of  leakage  through  these 
holes . 

2.2.8  Damage  to  Grout  Seals  and  Aquifer  Contamination 

The  results  from  the  Salmon  event  indicated  that,  as 

predicted,  neither  grout  seal  damage  nor  aquifer  contamina¬ 
tion  occurred  (Reference  2.11).  Because  it  was  predicted 
that  cracking  and  radioactive  injection  into  cracks  caused 
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3000* 


Figure  2.5  Sterling  Stemming  Plan 
(Redrawn  for  Reproduction  from  Reference  2.6) 
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by  the  Sterling  event  would  not  exceed  the  effects  caused  by 
Salmon  ,  no  hazards  to  grout  seals  or  aquifers  were  antici¬ 
pated  . 

2.2.9  Containment  Evaluation 

The  depth  of  burial  of  828  meters  provided  a  scaled 
depth  of  1,088  meters.  Peak  cavity  pressure  was  assumed  to 
be  less  than  overburden  pressure.  Even  the  largest  theoreti¬ 
cally  calculated  pressure  would  exceed  overburden  only 
slightly  and  for  only  a  few  milliseconds.  With  complete 
stemming,  therefore,  gross  dynamic  venting  appeared  unlikely. 

Minor  leakage  also  appeared  unlikely.  The  device  was 
under  366  meters  of  impermeable  salt  and  significant  fissur- 
ing  was  not  anticipated.  Illustrative  of  the  impermeability 
of  the  medium  was  the  vacuum  maintained  in  the  Salmon  cavity 
up  until  drill-back  several  months  after  detonation.  Ac¬ 
cordingly,  the  chance  of  escape  of  cavity  gases  through  the 
halite  was  considered  to  be  extremely  remote.  Stemming, 
which  was  continuous  to  the  surface  and  included  concrete 
plugs,  was  considered  to  be  adequate.  Assuming  effective 
grouting  of  the  casing,  the  possibility  of  leakage  appeared 
small . 

On  the  basis  of  the  foregoing  evaluations  and  the  as¬ 
sumption  that  construction  and  emplacement  would  be  accom¬ 
plished  in  accordance  with  the  proposed  plans,  it  was 
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predicted  that  the  Sterling  event  with  a  yield  of  44  kt 
would  be  contained. 

2.3  SEISMIC  PREDICTIONS 

2.3.1  Transmission  Model 

As  it  was  located  in  the  Salmon  post-shot  cavity,  the 
Sterling  event  was  expected  to  generate  seismic  waves  which 
would  be  transmitted  through  the  same  regional  geologic  en¬ 
vironment  as  those  generated  by  the  Salmon  event.  With  the 
exception  of  the  immediate  vicinity  of  the  source  where  a 
cavity  existed  and  the  surrounding  salt  had  undergone  some 
degree  of  alteration,  the  physical  properties  of  this  en¬ 
vironment  were  considered  to  be  the  same  as  at  the  time  of 
the  Salmon  event.  To  evaluate  fully  the  effect  of  these 
differences  would  have  required  a  detailed  knowledge  of  the 
degree  of  alteration  beyond  the  cavity  wall.  Even  if  this 
information  had  been  available,  the  effort  required  to  evalu¬ 
ate  the  effects  would  have  been  out  of  proportion  to  the 
total  effort  necessary  to  evaluate  Sterling  in  terms  of  safety 
hazards.  Therefore,  it  was  assumed  that  Salmon  experience 
was  directly  applicable  to  the  Sterling  event  except  for  the 
decoupling  factor.  The  Salmon  predictions  fitted  the 
observed  data  for  that  event  quite  well  (Reference  2.11). 
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The  Salmon  predictions  utilized  cube  root  scaling  relation¬ 
ships  applied  to  Gnome  subsurface  data  with  application  of 
transmission  coefficients  to  account  for  the  local  geology 
and  to  convert  subsurface  predictions  to  surface  predictions. 
Because  Salmon  surface  data  were  available,  cube  root  rela¬ 
tionships  were  used  to  predict  Sterling  surface  motions,  as¬ 
suming  fully  tamped  conditions. 

2.3.2  Prediction  Equations 

The  equations  discussed  below  were  based  on  a  fully 
tamped  situation.  To  predict  motions  for  the  Sterling  event 
under  decoupled  conditions  all  predictions  were  divided  by 
the  decoupling  factor  of  20.  This  factor  is  normally  applied 
to  peak  particle  displacement  but  was  utilized  here  to  pro¬ 
vide  predictions  for  all  motions,  i.e.,  peak  particle  sur¬ 
face  displacement,  velocity  and  acceleration. 

Acceleration 

The  cube  root  scaling  relationship  for  acceleration  is 


or  algebraically  rearranged 

a  =  K  W1/3(n_1)  R~n  (2) 

The  regression  equation  through  the  Salmon  peak  surface 
acceleration  data  was 

a  =  4.07  x  10^R“1#^5  (Reference  2.11)  (3) 
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Substituting  the  Salmon  yield  of  5.3  kt  into  Equation  (2) 
and  letting  4.07  x  10^  =  K  wl/3(n-l)  gives 

a  =  2.39  x  106  W-32  R-1'95  (4) 

which  was  the  equation  used  to  predict  surface  motions  from 
Sterling  in  a  fully  tamped  condition.  In  the  above  equations 
a  =  peak  surface  particle  acceleration  in  g 
W  =  yield  in  kilotons 
R  =  slant  distance  in  meters 
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which  was  used  to  predict  the  ground  motions  for  the  tamped 
case . 

Velocity 

Peak  particle  velocity  was  predicted  for  Salmon  using 
u  =  3.26  x  106  W55  R"1-64  (9) 

where 

u  =  peak  surface  particle  velocity  in  cm/sec 

This  equation  was  derived  from  Gnome  data  using  cube  root 
relationships  similar  to  those  just  discussed.  This  equation 
fitsthe  Salmon  data  sufficiently  well  that  no  modification  was 
required  for  predicting  the  Sterling  velocity  for  the  tamped 
situation . 

2.3.3  Predicted  Ground  Motions 

Ground  motion  predictions  were  made  only  for  the  maximum 
yield  of  0.44  kt  because  the  difference  in  motions  antici¬ 
pated  from  the  design  yield  of  0.36  kt  were  not  significantly 
less.  The  main  uncertainty  in  predicting  ground  motions  from 
the  Sterling  event  was  the  degree  of  decoupling  which  would 
occur.  The  Sterling  Technical  Concept  (Reference  1.2)  esti¬ 
mated  decoupling  might  vary  from  a  factor  of  20  to  a  factor 
of  160, depending  on  the  alteration  of  the  cavity  wall  which 
was  caused  by  the  Salmon  explosion.  The  only  previous 
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experiment  conducted  to  test  decoupling  was  Project  Cowboy, 
an  HE  experiment  at  Winnfield  Dome,  Louisiana.  The  Cowboy 
data  indicated  that  decoupling  was  achieved  with  HE  explosives. 
However,  the  theory  had  never  been  tested  with  nuclear  ex¬ 
plosives  . 

To  provide  ground  motion  predictions  which  would  be  use¬ 
ful  from  the  standpoint  of  safety,  the  motions  anticipated 
from  a  tamped  0.44  kt  explosion  were  calculated  using  Equations 
4.  8,  and  9  above.  The  predictions  for  selected  points  of 
interest  are  listed  in  Table  2.3.  These  predictions  were  made 
in  order  to  provide  an  upper  limit  on  the  ground  motions  from 
the  Sterling  event.  To  account  for  the  decoupling  effects  and 
to  provide  more  realistic  estimates  of  ground  motion,  the 
smallest  decoupling  factor,  20,  given  in  Reference  1.1  was 
used.  All  predictions  for  the  tamped  case  were  divided  by  20 
and  the  resulting  predictions  given  in  Table  2.4. 

Figures  2.6  through  2.8  show  peak  predicted  motions  versus 
slant  distance.  Each  graph  shows  the  predictions  for  Sterling 
assuming  a  decoupling  factor  of  20  and  assuming  a  fully  tamped 
(no  decoupling)  case. 

2.3.4  Predicted  Distances  to  Significant  Ground 
Motions 

Using  the  decoupling  factor  of  20,  it  was  predicted  that 
ground  motion  would  not  exceed  0.1  g  acceleration  beyond  1100 
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TABLE  2.3  PREDICTED  PEAK  SURFACE  MOTIONS,  0.44  kt  FULLY  TAMPED 


Slant 

Predicted  Motions 

Location 

Distance 
(meters ) 

Acceleration 
. . (9 ) 

Displacement 

(cm) 

Velocity 
( cm/ sec) 

1  mile  from  SZ 

1,810 

8. OxlO-1 

1 . 9xl0-1 

9 . 5x10° 

2  miles  from  SZ 

3,340 

2.4X10-1 

7 . 3xl0“2 

3.5x10° 

3  miles  from  SZ 

4,900 

1 . 2xl0-1 

4 . OxlO-2 

1.8x10° 

Baxterville 

6,000 

7.7xl0~2 

2. 9xl0-2 

1 . 3x10° 

Baxterville  Oil 
Field 

9,600 

3. lxlO-2 

1 . 4xl0-2 

5. 9x10-! 

Purvis 

14,500 

1.4xl0“2 

7 . 3xl0-3 

3. OxlO"1 

Lumberton 

17,700 

9.4xl0“3 

5.2xl0-3 

2. lxlO-1 

Columbia 

27,800 

3 . 9xl0-3 

2.6xl0“3 

1 . OxlO-1 

Hattiesburg 

32,000 

2. 9x10 “3 

2. OxlO"3 

8.0xl0-2 
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TABLE  2.4  PREDICTED  PEAK  SURFACE  MOTIONS,  0.44  kt  WITH  DECOUPLING  OF  20 


Slant  Predicted  Motions 


Location 

Distance 
(meters ) 

Acceleration 

(9) 

Displacement 

(cm) 

Velocity 
( cm/sec ) 

1  mile  from  SZ 

1,810 

4. 0xl0-2 

9.5xl0-3 

4 . 7xl0-1 

2  miles  from  SZ 

3,340 

1.2xl0-2 

3.6xl0-3 

1 . 8xl0-1 

3  miles  from  SZ 

4,900 

5 . 9xl0-3 

2. OxlO-3 

9. OxlO-2 

Baxterville 

6,000 

3.9xl0-3 

1 . 4x10 -3 

6 . 5xl0-2 

Baxterville  Oil 
Field 

9,600 

1 . 6xl0~3 

6 . 9xl0-4 

3. OxlO-2 

Purvis 

14,500 

7. OxlO-4 

3 . 6xl0~4 

1 . 5xl0-2 

Lumber ton 

17,700 

4. 7xl0-4 

2. 6xl0-4 

l.lxlO-2 

Columbia 

27 , 800 

2. OxlO-4 

1 . 3xl0-4 

5.0xl0-3 

Hattiesburg 

32,000 

1 . 5xl0-4 

1 . OxlO-4 

4. 0xl0-3 
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Slant  Distance,  meters 

Figure  2.7.  Predicted  Peak  Surface  Particle 
Displacement  versus  Slant  Distance,  0.44  kt. 


Slant  Distance,  meters 

Figure  2.8.  Predicted  Peak  Surface  Particle 
Velocity  versus  Slant  Distance,  0.44  kt. 


meters  slant  distance,  which  is  a  horizontal  distance  of  730 
meters  from  surface  zero.  If  no  decoupling  occurred,  0.1  g 
was  predicted  to  extend  to  a  slant  distance  5.3  km  (5.1  km 
from  surface  zero) . 

The  predicted  distance  to  a  particle  velocity  of  1  cm/sec 
was  1150  meters  slant  distance  or  800  meters  horizontal  from 
surface  zero  for  Sterling  using  a  decoupling  factor  of  20. 

If  no  decoupling  was  achieved  and  Sterling  effects  were  simi¬ 
lar  to  a  tamped  event,  the  distance  to  1  cm/sec  was  predicted 
to  be  7.0  kilometers. 

The  distance  to  0.001  g  for  the  Sterling  event  was  pre¬ 
dicted  to  be  12  km  assuming  a  decoupling  factor  of  20.  If 
no  decoupling  occurred,  peak  accelerations  of  0.001  g  were 
predicted  to  occur  as  far  as  55  km  from  the  Sterling  surface 
zero.  These  predicted  distances  were  average  values.  Azimuthal 
variation  in  ground  motion  was  observed  for  the  Salmon  event 
and  was  expected  to  occur  for  Sterling.  However,  it  was 
though l  that  the  variations  observed  at  Salmon  would  not  be 
applicable  to  Sterling  because  of  the  close-in  environmental 
differences.  Consequently,  no  attempt  was  made  to  incorpo¬ 
rate  the  observed  Salmon  azimuthal  variation  into  the  Sterling 
predictions . 
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The  pertinent  distances  are  summarized  below  and  shown 
graphically  with  respect  to  surrounding  towns  on  Figures  2.9 


and  2.10. 

Motion 
0.1  g 
1  cm/sec 
0.001  g 


Distance  from 

Surface  Zero,  km 

Fully  Tamped 

Decoupling  of  20 

5.1 

0.73 

7.0 

0.80 

55 

12 

2.3.5  Hazards  to  Subsurface  Facilities 

There  was  no  known  damage  to  any  of  the  oil  or  gas  wells 
in  the  vicinity  or  to  the  closest  pipeline  as  a  result  of  the 
Salmon  event.  Therefore,  it  was  concluded  that  Sterling,  at 
a  much  smaller  yield  and  decoupled,  would  not  present  any 
hazard  to  these  facilities. 


2 . 4  INSTRUMENTATION 

Seismic  instrumentation  was  installed  and  operated  by 
the  Special  projects  Party  of  the  U.S.  Coast  and  Geodetic 
Survey.  Instrumentation  consisted  of  1)  Accelerogr aphs  which 
record  three  orthogonal  components  of  particle  acceleration 
and  displacement  on  photographic  paper,  2)  National  Geophysi¬ 
cal  Co.  Type  21  seismometers  (NC-21)  which  record  three  com¬ 
ponents  of  particle  velocity  on  magnetic  tape  and, 3)  Wood-Anderson 
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Figure  2.10.  Predicted  Distance  to  Significant  Motions, 
Sterling  Event  Fully  Tamped 
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displacement  meters  which  record  two  horizontal  components 
of  particle  displacement  on  paper. 

The  accelerogr aphs  vrere  located  at  ten  nearby  structures 
within  a  radius  from  surface  zero  of  about  six  kilometers. 
Location  of  these  stations  are  shown  on  Figure  2.11. 

The  velocity  meters  and  Wood-Ander son  displacement  meters 
were  placed  at  selected  points  of  interest  such  as  nearby 
towns  and  cities.  The  locations  are  shown  on  Figure  2.12. 

At  two  stations,  10S  and  20S,  the  three-component  set 
of  NC -21  seismometers  were  supplemented  by  three  additional 
vertical  seismometers  at  1000  foot  intervals.  At  station  10S 
the  vertical  instruments  were  on  a  radial  line  in  a  direction 
awa y  from  ground  zero  while  at  Station  20S  the  vertical  in¬ 
struments  were  on  a  radial  line  toward  ground  zero.  Further 
details  concerning  the  instrumentation  for  the  Sterling  event 
may  be  found  in  Reference  2.12. 
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Figure  2.11  Map  Showing  Accelerograph  Station 
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Figure1  2.12  Area  Map  Showing  Additional  Station  Locations 
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CHAPTER  3 


RESULTS 

3 . 1  PHENOMENOLOGY  AND  CONTAINMENT 

The  Sterling  event  was  completely  contained.  Because 
there  has  been  no  post-shot  drilling  program  to  date,  no 
other  data  are  available. 

3.2  SEISMIC 

Seismic  signals  from  the  Sterling  event  were  so  weak 
that  many  of  the  instruments  failed  to  record  usable  data. 

The  best  data  are  those  recorded  on  magnetic  tape  with  the 
NC-21  velocity  meters.  Of  the  eight  NC-21  stations  in  opera¬ 
tion,  seven  recorded  usable  data,  although  in  some  cases  not 
all  components  (at  a  station)  were  usable.  No  usable  data 
were  obtained  at  Hattiesburg. 

The  NC-21  velocity  meter  data  were  processed  by  analog 
computer.  Corrections  were  made  for  the  frequency  response 
characteristics  of  the  instruments  to  obtain  corrected 
velocity.  Acceleration  and  displacement  traces  were  de¬ 
rived  by  integration  and  differentiation  of  the  corrected 
velocity  traces.  The  velocity  traces  were  also  filtered 
with  narrow  band-pass  filters  at  pre-selected  center  frequen¬ 
cies.  The  peak  velocities  at  each  frequency  were  plotted 
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versus  frequency  to  provide  amplitude-frequency  relation¬ 
ships  at  each  station.  The  center  frequencies  were  selected 
so  as  to  include,  where  possible,  all  frequencies  which 
were  of  significant  amplitude.  At  three  stations;  Columbia, 
Lumberton,  and  Bass  Memorial  School  at  Purvis;  high  cutoff 
filters  of  16  hz  were  used  during  the  field  recording  and 
at  both  Baxtervi'lle  stations  the  cutoff  was  at  23  hz. 

Band-pass  filtering  was  not  carried  beyond  these  frequen¬ 
cies.  In  general,  the  peak  amplitude  at  these  stations 
appeared  to  be  occurring  at  frequencies  lower  than  the  high 
cutoff  field  filters,  so  it  was  possible  to  identify  the 
significant  frequencies . 

Of  the  ten  Accelerograph  stations,  only  four  produced 
traces  with  sufficient  amplitudes  for  analysis.  These 
were  recorded  with  the  accelerometers.  The  Carder  dis¬ 
placement  meters  (the  Accelerograph  consists  of  three 
accelerometers  and  three  Carder  displacement  meters)  did 
not  record  any  motion.  The  useable  accelerometer  traces 
exhibited  very  high  frequency  motions,  generally  well 
outside  the  flat  portion  of  the  response  curve.  These 
traces  were  digitized  and  corrected  for  the  frequency 
response  characteristics  of  the  instruments.  Particle 
velocity  traces  were  derived  from  these  corrected  acceleration 
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data  and  were  then  band-pass  filtered  by  digital  methods 
to  provide  amplitude-frequency  information  comparable  to 
that  obtained  from  the  NC-21  data.  No  readable  motions 
were  obtained  at  the  Wood-Anderson  displacement  meter 
stations  0 

The  peak  values  of  ground  motion  are  listed  in  Tables 
3.1  through  3.3C  The  tables  include  the  peak  value  for 
each  component  of  motion  and  the  peak  resultant  vector 
value o  The  resultant  vectors  were  obtained  by  analyzing 
simultaneously  the  three  components  of  motion  at  a  station 
in  order  to  determine  the  absolute  value  of  ground  motion. 
The  resultant  vectors  are  the  instantaneous  value  of  the 
square  root  of  the  sum  of  the  squares  of  the  individual 
components.  Where  only  two  components  of  motion  were  re¬ 
corded,  peak  resultant  vectors  were  also  calculated.  The 
peak  resultant  vector  values  are  plotted  versus  slant 
distance  on  Figures  3.1  through  3*3*  Least  squares  re¬ 
gression  lines  were  fitted  to  the  data  and  the  standard 
errors  of  estimate  calculated*,  These  are  given  on  the 
figures  0 

The  amplitude-frequency  curves,  derived  by  filtering 
the  velocity  signals,  are  shown  on  Figures  3*4  through 
3.16. 
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TABLE  3.1.  PEAK  SURFACE  ACCELERATION 


Slant 
Distance 
( meters ) 

Peak  Acceleration,  g 

Station 

Vertical 

Radial 

Transverse 

Resultant 

1 

2,550 

3.18xl0~2 

1.35x10“2(N-S) 

3 . 18x10” 2 

2 

2,620 

2„20xl0“2 

9.25x10“3(N-S) 

1 . 68xl0“2(E-W) 

2 . 28xl0“2 

3 

5,970 

— 

— 

— 

_ 

4 

1,830 

8. 6 9x10“ 2 

3.10x10“2(N-S) 

3 . 03x10“ 2 (E-W ) 

8.99xl0“2 

5 

5,210 

— 

— 

_ 

_ 

6 

4,720 

— 

— 

_ _ 

_ 

7 

4,320 

— 

— 

_ 

—  - 

8 

2,060 

2. 38x10“ 1 

8.81x10“2(N-S) 

— . 

2. 38x1 0"1 

9 

5,390 

— 

- - 

_ 

10 

3,520 

— 

— 

_ 

. 

Baxter ville 

Post  Office* 

6,100 

1.61xl0"3 

7 . 95xl0“4 

6 . 35xl0“4 

1 . 69xl0”3 

School* 

6,500 

1 . 43xl0“3 

8 . 74xl0-4 

5 . 13xl0“4 

1 . 55xl0“3 

Purvis  (Bass 

16,000 

1 . 55xl0“4 

4. 94x1 0“5 

4.59xl0“5 

1 ,56xl0“4 

Memorial 

School)* 

10  S* 

17,800 

2.22x10 

A 

7.28xl0“b 

3.97xl0“5 

2 . 22xl0“4 

1.70xl0“4 
2. 02x10“ 4 
1.50xl0“4 


Lumberton* 

19,000 

3.85xl0“5 

4 . 65xl0“5 

2.69xl0“5 

5. 15xl0“5 

Columbia* 

26,000 

— 

— 

_ 

20  S* 

31,700 

1 . 56xl0“5 

1 .72xl0“5 

2 . 31xl0“5 
1.84xl0“5 

1 .97xl0“5 

1.06xl0“5 

2. 12xl0“5 

Hattiesburg* 

32,700 

1 

*  Derived  from  corrected  velocity. 
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TABLE  3.2.  PEAK  SURFACE  DISPLACEMENT 


Slant 

Peak  Displacement,  cm 

Station 

Distance 
( meters ) 

Vertical 

Radi al 

Transverse 

Resultant 

1 

2,550 

, 

2 

2,620 

— 

— 

— 

— 

3 

5,970 

— 

— 

__ 

— 

4 

1,830 

— 

— 

— 

— 

5 

5,210 

* — 

— 

— 

— 

6 

4,720 

— 

— 

— 

— 

7 

4,320 

— 

— 

— 

— 

8 

2,060 

— 

— 

— 

— 

9 

5,390 

— 

— 

— 

— 

10 

3,520 

— 

— 

— 

— 

Baxter ville 

Post  Office* 

6,100 

2. 57x1 O-4 

2.26X10"4 

_ 

3.08X10"4 

School* 

6,500 

3 . 23xl0-4 

1 . 83xl0-4 

7 . 16x10“"’ 

3 . 24xl0~4 

Purvis  (Bass 

16,000 

8 . 7  8xl0-5 

3.45X10"5 

3.45xl0-5 

1 .02xl0"4 

Memorial 

School ) * 

10  S* 

17 , 800 

8 . 50x10"-’ 

9.24xl0-5 

2 . 49x10“'’ 

9 . 57xl0~4 

Lumber ton* 

19,000 

8.08X10"5 

1 . 13xl0~4 
7.23xl0~5 

4 . 90x10"'’ 

6 . 21xl0-5 

3.09xl0-5 

7.11xl0-5 

Columbia* 

26,000 

— 

— 

— 

— 

20  S* 

31,700 

— 

— 

— 

— 

Hattiesburg* 

32,700 

-* 


Derived  from  corrected  velocity. 


TABLE  3.3.  PEAK  SURFACE  VELOCITY 


Slant 
Distance 
(meters ) 

Peak  Velocity,  cm/sec 

Station 

Vertical 

Radial 

Transverse 

Resultant 

1* 

2,550 

1. 65x10" j 

8.06xl0"2(N-S) 

1. 65x10" 1 

2* 

2,620 

1. 27x10" 1 

5.95x10  (N-S) 

8.05x10“2(E-W) 

1. 34x10" 1 

3 

5,970 

— 

— 

— 

-  ■ 

4* 

1,830 

3. 48x10" 1 

1. 37x10" 1( N-S) 

1. 27x10" 1(E-W) 

3.58X10"1 

5 

5,210 

— 

— 

_ 

—  - 

6 

4,720 

— 

— 

_ 

,  -  ■ 

7 

4,320 

— 

- - 

8* 

2,060 

6 . 09x1 0-1 

1. 50x10" X( N-S) 

— 

6. 09x1 0"1 

9 

5,390 

— 

_ 

■ 

10 

3,520 

— 

— 

_ 

Baxterville 
Post  Office 

6,100 

1.59xl0"2 

1.24xl0"2 

7. 06x10“ 3 

1.6 8x10" 2 

School 

6,500 

1.64x10““ 

1 . 31xl0"2 

5.34xl0"3 

1.69xl0"2 

Purvis  (Bass 

16,000 

3. 08x10" 3 

8.97xl0“4 

7.59xl0"4 

3 . 15x1 0“3 

Memorial 
School ) 

10S 

17,800 

2.81xl0“3 

1.62xl0"3 

5. 98x10“ 4 

2 . 85x1 0"3 

2  .  07x10" 3 
2 .41x10" 3 
1.74xl0"3 


Lumberton 

19,000 

1.00x10" 3 

1 . 35xl0“3 

7. 59x10" 4 

1.59x10" 3 

Columbia 

26 , 000 

— 

— 

4.03xl0“4(E-W) 

- 

2  OS 

31,700 

5. 48x1 0"4 
4. 90x1 0"4 

5 . 52xl0"4 
5. 90x10" 4 

6.83xl0“4 

2 . 40xl0"4 

6.93xl0"4 

Hattiesburg 

32,700 

*  Derived  from  corrected  acceleration. 
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Peak  Acceleration 
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Figure  3.1  Observed  Peak  Resultant  Vectors  of 
Surface  Par  tide  Accelcr  at  ion 
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Figure  3.2  Observed  Peak  Resultant  Vector 
of  Surface  Particle  Displacement 


Slant  Distance,  meters 


Figure  3.3  Observed  Peak  Resultant  Vectors 
of  Surface  Particle  Velocity 


Peak  Particle  Velocity,  cm/sec 


Frequency,  hz 

Figure  3.4  Peak  Rand-Pass  Filtered  Particle  Velocity 
versus  Frequency,  Station  4 
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Frequency,  hz 


Figure  3.5  Peak  Band-Pass  Filtered  Particle  Velocity 
versus  Frequency,  Station  8 


Frequency,  hz 


Figure  3.6  Peak  Band-Pass  Filtered  Particle  Velocity 
versus  Frequency,  Station  1 


Peak  Particle  Velocity,  cm/sec 


Figure  3.7  Peak  Band-Pass  Filtered  Particle  Velocity 
versus  Frequency,  Station  2 
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Peak  Particle  Velocity,  cm/sec 


Peak  Particle  Velocity,  cm/sec 


Frequency,  hz 

Figure  3.14  Peak  Band-Pass  Filtered  Particle  Velocity 

versus  Frequency,  Columbia 


Peak  Particle  Velocity,  cm/sec 


Frequency,  hz 

Figure  3.15  Peak  Band-Pass  Filtered  Particle  Velocity 
versus  Frequency,  Station  20S 
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CHAPTER  4 


ANALYSIS  AND  DISCUSSION 

4.1  PHENOMENOLOGY  AND  CONTAINMENT 

As  stated  in  Chapter  3,  the  only  information  available 
is  that  the  event  was  completely  contained  as  predicted. 

The  prediction  that  containment  would  be  achieved  was  based 
on  the  prediction  and  evaluation  of  the  individual  effects 
and  their  interrelation  in  terms  of  containment.  Because 
Sterling  was  contained,  there  is  a  qualitative  substantia¬ 
tion  of  ERC's  evaluation  of  the  adequacy  of  stemming  and 
the  depth  of  burial.  However,  none  of  the  other  predicted 
quantities  can  be  compared  with  results,  because  there  has 
been  no  post-shot  exploratory  program  to  date. 

4.2  SEISMIC 

Comparisons  of  predictions  and  results  are  shown  on 
Figures  4.1  through  4.3  for  acceleration,  displacement,  and 
velocity,  respectively.  Each  figure  shows  the  predictions, 
for  0.35  kt,  for  a  fully  tamped  situation  and  for  a  decoupling 
factor  of  20.  For  comparison  with  the  predictions,  the 
least  squares  regression  lines  which  were  fitted  to  the  peak 
resultant  vectors  of  motions  are  shown.  In  general,  the 
observed  motions  from  the  Sterling  event  were  associated 
with  high  frequencies,  and  the  peak  amplitudes  attenuated 
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with  distance  much  more  rapidly  than  did  the  peak  motions 
observed  for  the  Salmon  event. 

4.2.1  Acceleration 

The  comparison  between  predicted  and  observed  ac¬ 
celeration  is  shown  on  Figure  4.1.  The  regression  equation 
through  the  data  shows  an  attenuation  rate  proportional  to 
£-3.04^  an  extremely  rapid  attenuation  when  compared  with 
the  predicted  rate  of  R"1*95  which  was  based  on  Salmon  data. 
The  predictions  using  a  decoupling  factor  of  20  are  lower 
than  the  observed  data  at  the  close  distances  and  higher 
at  the  long  distances.  The  scatter  in  the  observed  data 
was  measured  by  calculating  the  standard  error  of  estimate, 
a.  For  the  acceleration  data,  the  value  of  a  was  1.79 
which  indicates  that  68%  of  the  data  points  are  within  a 
factor  of  1.79  of  the  regression  equation,  assuming  nor¬ 
mal  distribution.  By  comparison,  the  standard  error  cal¬ 
culated  for  the  Salmon  acceleration  data  was  1.66  (Reference 
2.11) . 

4.2.2  Displacement 

The  displacement  comparison  is  shown  on  Figure  4.2. 

The  predictions,  using  a  decoupling  factor  of  20,  were 
higher  than  observed  by  factors  ranging  from  about  4  at 
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Acceleration 


placement 


the  closest  station  to  about  2,5  at  the  most  distant.  The 
data  showed  very  little  scatter;  the  standard  error  of  esti¬ 
mate  was  a  factor  of  1.11  from  the  regression  line.  The 
attenuation  rate  indicated  by  the  regression  line  through 
the  data  was  R"-1-*25  compared  with  the  predicted  rate  of 
p-1.60.  However,  only  five  displacement  vectors  were  obtained. 
Therefore,  the  indicated  scatter  and  the  rate  of  attenuation 
of  peak  displacement,  calculated  statistically  on  a  very 
small  sample,  are  not  necessarily  representative  of  dis¬ 
placement  behavior  in  general.  For  example,  the  Salmon 
data  points  showed  a  standard  error  of  1.66  (Reference  2.11). 

4.2.3  Velocity 

The  comparison  between  predicted  and  observed  particle 
velocity  is  shown  on  Figure  4.3.  The  scatter  in  the  velocity 
data  was  small.  The  standard  error  was  a  factor  of  1.40 
compared  with  the  factor  of  1.66  associated  with  the  Salmon 
data.  The  predictions  using  the  decoupling  factor  of  20 
were  higher  than  the  observed  velocities  by  factors  ranging 
from  about  1.1  to  6,  from  the  closest  to  the  most  distant 
stations,  respectively.  The  observed  data  showed  a  rate  of 
attenuation  with  distance  of  R-2*27  compared  with  the  pre¬ 
dicted  rate  of  R 
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4.2.4  Azimuthal  Distribution  of  Peak  Amplitudes 


The  ground  motions  recorded  from  the  Salmon  event  ex¬ 
hibited  an  asymmetrical  distribution  with  azimuth  (Reference 
2.11).  For  comparison,  the  Sterling  data  were  plotted  and 
contours  were  drawn  showing  the  azimuthal  distribution  of 
the  peak  amplitudes  measured  from  the  Sterling  event.  Peak 
vertical,  radial  and  transverse  components  of  particle 
velocity  are  shown  on  Figures  4.4  through  4.6  Peak  ac¬ 
celerations  are  shown  on  Figures  4.7  through  4.9. 

The  azimuthal  variations  observed  in  the  Sterling  data 
were  relatively  minor  compared  with  Salmon.  The  high  motions 
at  Station  8,  about  2  km  west  of  surface  zero,  indicated 
preferential  propagation  to  the  west.  However,  this  might 
have  been  a  local  phenomenon  because  the  one  component  of 
motion  (transverse  velocity)  measured  at  Columbia  (see  Fig¬ 
ure  4.6)  appeared  normal  with  respect  to  the  motions  at 
other  stations.  The  distribution  of  the  peak  velocities 
and  accelerations  was  generally  similar  if  the  same  components 
were  compared,  but  differences  were  noted  between  the  dis¬ 
tribution  of  vertical  and  horizontal  components.  The 
vertical  components  of  both  velocity  and  acceleration  were 
slightly  low  at  Lumberton  indicating  increased  attenuation 
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Figure  4.7  Azimuthal  Distribution  of  Peak  Vertical  Particle  Acceleration 
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of  the  vertical  motion  in  the  southeast  direction.  How¬ 


ever,  the  horizontal  components  of  velocity  and  accelera¬ 
tion  were  about  normal  at  Lumberton  but  low  at  Purvis  to 
the  east,  and  Baxterville  and  10S  to  the  south. 

4.2.5  Amp 1 i tude -Frequency 

The  amplitude-frequency  curves  shown  on  Figures  3.4 
through  3.15  indicate  that  the  peak  amplitudes  of  ground 
motion  were  associated  with  relatively  high  frequencies. 

The  closest  stations  (Figures  3.4  through  3.7)  show 
that  the  peak  velocities  were  associated  with  very  high 
frequency  energy.  The  amplitude-frequency  curves  for 
Stations  1,  2  and  4  reach  a  maximum  between  20  and  40  hz , 
while  at  Station  8  the  peak  occurs  at  about  80  hz . 

With  increasing  distance  the  frequency  at  which  the 
peak  amplitudes  occur  tends  to  decrease.  For  example 
Baxterville  Post  Office  (6.1  km),  Baxterville  School 
(6.5  km),  Bass  Memorial  School  near  Purvis  (16  km),  and 
10S  (17.8  km)  show  peak  amplitudes  associated  with  fre¬ 
quencies  between  10  and  20  hz . 

Lumberton  and  Columbia,  19  and  26  km,  respectively, 
received  peak  motions  with  frequencies  between  3.5  and  10 
hz  while  at  Station  20S,  31.7  km,  the  peak  motions  were 
at  frequencies  between  3  and  5  hz . 
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The  decrease  in  frequency  with  increasing  distance  is 
illustrated  on  Figure  4.10.  The  frequency  of  the  largest 
band-pass  filtered  velocity  is  shown  versus  distance  for 
each  component.  In  a  few  cases,  two  frequencies  had  equal 
amplitudes  so  the  average  frequency  was  plotted.  Least- 
squares  regressions  through  these  data  show  that  the 
change  in  frequency  was  proportional  to  R“*^  on  the 
vertical  component,  R  on  the  radial,  and  R  on 

the  transverse. 

4.2.6  Observed  Decoupling  Factors 

The  level  of  the  Sterling  ground  motions  and  the  high 
rates  of  attenuation  for  particle  acceleration  and  velocity 
(Figures  4.1  and  4.3)  considered  together  with  the  frequency- 
distance  relationships  (Figure  4.10)  indicate  that  the 
Sterling  seismic  signal  was  decoupled  and  that  decoupling 
was  frequency-dependent.  A  direct  measure  of  the  decoupling 
as  a  function  of  frequency  results  from  comparison  of  the 
Salmon  and  Sterling  band-pass-filtered  particle  velocities 
at  Stations  10S  and  20S  where  identical  instrumentation  was 
employed  for  both  events.  Five  traces  are  available  for 
comparison  at  Station  10S ;  one  radial,  one  transverse  and 
three  vertical.  Station  20S  provided  a  radial,  a  trans¬ 
verse,  and  two  vertical  traces. 
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Figure  4.10  Frequency  of  the  Peak  Band-Pass  Filtered 
Velocity  versus  Slant  Distance 

-76- 


The  Salmon  band -pass ~f iltered  velocities  at  each 
frequency  were  scaled  to  the  Sterling  yield  of  0.35  kt  using 
the  frequency-dependent  yield  scaling  exponents  from  the 
extrapolation  method  discussed  in  Reference  4.1.  These 
scaled  velocities  were  then  compared  at  each  frequency  with 
the  Sterling  velocities.  After  scaling  the  Salmon  data  to 
the  Sterling  yield,  the  Salmon/Ster ling  ratios  are  the  de¬ 
coupling  factors  for  particle  velocity.  Decoupling  factors 
were  calculated  for  fifteen  different  frequencies  ranging 
from  .41  to  29  hz .  These  values  are  listed  in  Table  4.1. 
Although  the  decoupling  factors  exhibit  some  scatter,  there 
are  no  noticeable  major  differences  between  the  factors  ob¬ 
served  at  the  two  stations  and  the  factors  were  averaged 
for  each  frequency.  Separate  averages  were  compiled  for  the 
vertical  and  horizontal  components.  These  average  factors 
are  plotted  versus  frequency  on  Figure  4.11.  The  maximum 
decoupling  (a  factor  of  about  200)  occurred  for  frequencies 
between  2  and  5  hz .  The  amount  of  decoupling  decreased 
exponentially  for  both  higher  and  lower  frequencies,  dropping 
to  a  factor  of  about  50  at  .41  hz  and  to  a  factor  of  about 
15  at  29  hz .  The  decoupling  exceeded  a  factor  of  100  for 
frequencies  between  1  and  8  hz . 

The  observed  dependency  of  decoupling  on  frequency 
provides  an  explanation  for  the  high  rates  of  attenuation 


TABLE  4.1.  OBSERVED  DECOUPLING  FACTORS  AT  STATIONS  10S  AND  2 OS 


f 

(hz) 

Vertical  Components 

Average 

Station  10S 

Station  20S 

Zo 

Z1000’ 

Z2000' 

Z1000' 

Z3000' 

.41 

48 

189 

57 

52 

57 

81 

.55 

30 

117 

44 

83 

81 

71 

.75 

80 

86 

70 

71 

126 

87 

1.0 

150 

99 

113 

151 

235 

150 

1.35 

125 

91 

120 

90 

223 

130 

1.85 

122 

144 

118 

141 

215 

148 

2.5 

269 

176 

177 

282 

358 

252 

3.4 

232 

224 

175 

276 

268 

235 

4.6 

217 

240 

175 

311 

191 

227 

6.2 

210 

188 

158 

310 

163 

206 

00 

57 

90 

63 

180 

88 

96 

11.3 

37 

66 

35 

58 

63 

52 

15.5 

17 

22 

12 

30 

22 

21 

21 

16 

29 

11 

26 

19 

20 

29 

7.4 

24 

9.4 

— 

14 

_ Horizontal  Components _ 

_ Station  IQS _  _ Station  2 OS 

Radial _ Transverse _ Radial  Transverse  Average 


.41 

24 

73 

.55 

36 

93 

.75 

52 

82 

1.0 

124 

149 

1.35 

127 

181 

1.85 

161 

163 

2.5 

228 

251 

3.4 

177 

266 

4.6 

112 

232 

6.2 

116 

190 

8.4 

105 

143 

11.3 

63 

66 

15.5 

25 

12 

21 

19 

14 

29 

13. 

16 

50 

20 

42 

51 

65 

61 

63 

35 

58 

62 

72 

102 

117 

118 

136 

215 

140 

170 

378 

282 

285 

255 

227 

231 

208 

149 

175 

165 

170 

160 

89 

130 

117 

79 

50 

65 

39 

55 

33 

39 

33 

26 

— 

15 
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Frequency,  hz 

Particle  Velocity  Decoupling  Factors  Calculated 
from  observed  Salmon  and  Sterling  Band  Pass 
Filtered  Velocities  at  Stations  10  S  and  20  S 


which  were  observed  for  the  Sterling  peak  motions.  On 
Figure  4.3,  the  peak  velocities  at  the  closest  distances 
are  below  the  fully  coupled  prediction  by  a  factor  of  about 
20.  These  stations  recorded  relatively  more  high  frequency 
energy  (30  to  50  hz)  than  more  distant  stations. 

As  distance  increases,  peak  amplitudes  of  particle 
velocity  occur  at  progressively  lower  frequencies.  More¬ 
over,  departures  of  recorded  peak  values  from  the  Salmon- 
derived  attenuation  rates  follow  a  similar  trend.  The  re¬ 
sults  of  this  analysis  (Table  4.1,  Figure  4.11)  show  the 
variation  in  decoupling  factor  with  frequency  of  the  peak 
amplitudes . 

In  summary,  the  particle  velocities  from  the  Sterling 
event  were  decoupled  by  as  much  as  a  factor  of  200,  but 
only  for  frequencies  in  the  2  to  5  hz  range.  The  observed 
decoupling  steadily  decreased  at  frequencies  outside  this 
range.  The  decoupling  factor  was  over  100  throughout  the 
1  to  8  hz  interval.  (See  Figure  4.11.) 

4.2.7  Comparison  of  Theoretical  and  Observed 
Decoupling  Factors* 

The  frequency  dependence  of  the  decoupling  factor  noted 
in  the  previous  section  suggested  a  study  of  the  decoupling 


*  This  section  contributed  by  J.  R.  Murphy 
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"theory  presented  by  Latter  et .  al .  ,  to  determine  if  the 
theory  is  compatible  with  the  observations  at  the  distances 
under  consideration. 

The  Fourier  transform  of  the  displacement  at  a  distance 
R  from  a  tamped  shot  is  given  by  (Reference  4.1): 


where  c  is  the  compressional  wave  velocity  in  the  medium, 
a'  is  the  radius  at  which  the  medium  begins  to  behave 
elastically,  p(u>)  is  the  Fourier  transform  of  the  pressure 
that  acts  at  a'  and  co^  =  c/a' .  With  R  »  c/co,  Equation  (1) 
reduces  to: 


rT  (w)  = 


ID '2  +  i(U  '  W  -  Al2 CD2 
O  O  4u> 


Similarly,  the  Fourier  transform  of  the  elastic  dis¬ 
placement  produced  by  a  step  pressure  p,  at  a  distance  R 
for  a  decoupled  shot  is  given  by  (Reference  4.2): 


((B)  =  -  Eg— 

v  8tt|J.R 


+1(0  QC0 


^+2M-  ^2 


where  a  is  the  radius  of  the  spherical  cavity  sufficiently 
large  that  the  pressure  developed  on  the  wall  does  not 
exceed  the  elastic  limit  of  the  medium  and  <0o  =  c/a. 
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Thus,  the  decoupling  factor  as  a  function  of  frequency 
is  given  by 


A 

h  (“) 


(®) 


We  now  calculate  this  decoupling  factor  for  several  physically 
realizable  pressure  profiles  at  R  =  a*  . 

Step  Function 

First,  let  the  pressure  profile  at  R  =  a'  be  given  by 
p  H  ( t )  where  H(t)  is  the  unit  step  function,  Then  p(to)  is 
given  by : 


0)  /  0,  p(w>)  -  . 

2ttcd  . 


Letting  X  =  u  and  substituting  for  p(ou)  in  Equation  (la); 

A  P  0 

^tC^)  =  aT  (o)q  +cuoia)o(J0  ~3/4fu2) 

f  (cu)  a(U)o2+iajotu  -3/4cu2) 

Setting  a  =  (mQ2  -  3/4u)2),  p  =  (w'0Z  -  3/4co2),  this  gives 
a  decoupling  factor 

aT  /cv2  +  tuo20)2 
a  /p2+cu o  2  a)2 


(3) 


A 

A  /  X 
!d(») 
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Interpolation  Function 

Next,  let  the  pressure  profile  at  R  =  a1  be  given  by 
p(t)  =  p  sine  t  H(t)  | 


p  sine  t  H(t) 


where 


sine  t  =  sxn.  ■TT.,.T. 


t  =  T 


t  =  —  .  Then  p(u>)  =  °  f 

To  2tt  J 


sine  t  H ( t ) e 


-icoT0T 


pTo  ni-g 


-±  leg  JZlH 

2"  ^  .  1/2 

2rr 


where  II  is  the  rectangle  function  given  by 

I0)T 

0,  -SL  >  1/2 

1,  I  !!s  I  <1/2 


Letting  n 

n  =  - 


*T0 

2rr 

+  1/2 

COT 

o 

2n 

-  1/2 

p(to)  =  1—2  ^ n  -  iy)  . 

2n 


Therefore, 

A 


fT(u))  _  al  To03  [(Y^-ncQQCQ^  +i(  ncY+vu)0(u )  ] 

f ‘  (a>)  a  (P+itu>) 
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This  gives  a  decoupling  factor 


(4) 


A 

fT(u>) 


?D<“> 


a'  T0cu  1/7YQ'-n(Ju0cu)2+(nQ'+YtJu0cu)2 


a  y/ P2-i-(uJL  2  (JO2 


Exponential  Function 

Finally,  let  the  pressure  profile  at  R  =  a’  be  given 


Then 

CO 

P(»)  =  —  f  e  ,T*  H( T )  e^icuT°T  dT 
2tt  _< 


_  PTp 
2tt 


1  -  io>T0 
1  +  (^T0)2 


Letting 


nT  = 


0)To 

Y*  ° 


1  +  (0>To)2  ,  l+(e«T0)2 


p(tu)  =  El£  (n»  -  i  Y») 
2tt 
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Therefore,  the  decoupling  factor  in  this  case  is  given 
by : 

ar  T0uu  /( y  T  cv  -  n  T  to  o0) ) 2  +  ( n  1  at  +y  T  t»0uT)  ^ 
a  /&2+w'0  2  u)2 

In  evaluating  Equations  (3)  through  (5),  the  following 
constants  were  used;  c  =  4267  m/sec,  a  =  16.8  m,  Tq  =  0.02 
seconds.  The  elastic  radius  of  Salmon,  af  ,  was  taken  to  be 
300  m  and  400  m  respectively  in  accordance  with  the  lower 
and  upper  bounds  determined  by  other  investigators  (Refer¬ 
ence  4.3). 

The  decoupling  factors  corresponding  to  Equations  (3), 

(4) ,  and  (5)  with  a!  =  300  m  are  shown  in  Figures  4.12,  4.13, 
and  4.14  respectively,  and  those  corresponding  to  the  same 
equations  with  af  =  400m  are  shown  in  Figures  4.15-4.17.  No  at¬ 
tempts  have  been  made  to  obtain  realistic  magnitudes  of  the  pres¬ 
sures  acting  at  a  and  a?,  and  therefore  the  magnitudes  appearing 
on  the  ordinates  are  not  necessarily  representative  of 

the  values  to  be  expected.  It  can  be  seen,  however,  that 
the  shapes  of  the  curves  corresponding  to  Equations  (4)  and 

(5)  are  very  similar  to  the  shape  of  the  decoupling  curve 
from  the  observed  BPF  data.  In  all  cases,  the  calculated 
curves  roll  off  at  the  high  frequency  end  as  does  the  observed. 


^(o)) 

A 

rD(w) 
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In  conclusion,  it  may  be  said,  that  the  theory  predicts 


a  frequency  dependence  of  the  decoupling  factor  which 
in  good  agreement  with  the  measured  data. 


is 
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Decoupling  Factor 
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Figure  4.13  Decoupling  Factor  for  Interpolation  Function 
Pressure  Profile,  p(t)  =  p  sincrH(r) 
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Figure  4.15  Decoupling  Factor  for  Step  Function  Pressure 
Profile  p(t)  =  p  H(t). 
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Figure  4.17  Decoupling  Factor  for  Exponential  Function 
Pressure  Profile,  p(t)=p  H(r). 


CHAPTER  5 


CONCLUSIONS 

5 . 1  CONTAINMENT 

The  prediction  that  the  Sterling  event  would  be  con¬ 
tained  was  substantiated.  Although  no  post-shot  drilling 
program  has  yet  been  conducted,  the  tact  that  Sterling  was 
contained  substantiated  ERC*s  evaluation  of  such  items  as 
the  adequacy  of  stemming  and  the  depth  of  burial. 

5.2  SEISMIC 

All  observed  peak  motions  were  less  than  predicted 
for  a  fully  tamped  event  of  this  yield.  Compared  with  the 
predictions  using  a  decoupling  factor  of  20,  the  observed 
peak  displacements  and  velocities  were  less  than  predicted 
at  all  stations,  and  peak  accelerations  were  larger  than 
predicted  at  stations  within  4  km  and  less  at  more  distant 
stations.  The  observed  rates  of  attenuation  of  particle 
acceleration  and  velocity  were  considerably  greater  than 
the  predicted  rates  which  were  based  on  Salmon  data. 

The  data  did  not  show  abnormal  scatter.  Some  azi¬ 
muthal  variations  were  observed,  but  the  variations  were 
small  in  comparison  to  the  observed  azimuthal  variations  in 
the  Salmon  data. 

The  frequency  content  was  defined  at  eleven  locations. 
Peak  amplitudes  were  associated  with  frequencies  ranging 


-93- 


from  3  to  80  hz.  A  general  decrease  in  frequency  with 
increasing  distance  was  observed. 

The  seismic  signal  from  the  Sterling  event  was  de¬ 
coupled  and  the  amount  of  decoupling  was  dependent  on 
frequency.  Between  2  and  5  hz,  the  particle  velocities  were 
decoupled  by  a  factor  of  about  200.  Between  1  and  8  hz, 
the  decoupling  factor  was  100  or  greater .  An  expontential 
decrease  in  decoupling  was  observed  for  frequencies  outside 
the  2  to  5  hz  range.  The  frequency  dependence  of  the  ob¬ 
served  decoupling  factors,  i.e.,  the  shape  of  the  decoupling 
factor  versus  frequency  curve,  agreed  very  closely  with 
calculations  made  using  the  decoupling  theory. 
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Table  A-l . 


Internal  Energy  Table  (Reference  2.9) 


T°K  E  k  cals/kg 


300 

51.18 

700 

122.6 

2,400 

513.6 

3,800 

1,083. 

4,200 

1,320. 

5 , 000 

1,772. 

6,000 

2,278. 

7,000 

3,025. 

8,500 

5,400. 

10,000 

7,699. 

11,000 

9,071. 

12,000 

10,070. 

14,000 

11,900. 

15,849 

13,420. 

17,783 

15,750. 

25,119 

28,500. 

28,184 

33,820. 

31,623 

29,080. 

35,481 

44,010. 

39,811 

50,260. 

44,668 

59,470. 

50,000 

72,000. 

70,000 

123,000. 

90,000 

186,000. 
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Figure  A1  (Redrawn  for  Reproduction  from  Fig.  D-2! 
Reference  2.9  with  permission  of  LRL) 
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List  of  Abbreviations  for  Technical  Agencies 
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